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r-?', ABSTRACT 

o I Recently there has been a hint of a gamma-ray line at 130 GeV originated from the galactic 
centre after the analysis of the Fermi-LAT satellite data. Being monochromatic in nature, it 
rules out the possibility of having its astrophysical origin, and there has been a speculation that 
this line could be originated from dark matter annihilation. In this work, we propose a two 
component dark matter scenario where an extension of the Standard model by an inert Higgs 

\Q I doublet and a gauge singlet scalar concocted with Z 2 x Z' 2 symmetry, is considered. We find that 
our scenario can not only explain the 130 GeV gamma-ray line through dark matter annihilation 
but also produce the correct dark matter relic density. We have used the Standard Model Higgs 

^ ■ mass around 125 GeV as intimated by the LHC data. 
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1 Introduction 



Although the existence of dark matter (DM) is now established by various astronomical mea- 
surements and observations where the gravitational effects of this huge amount of dark matter is 
manifested, the particle nature of the dark matter still remains unknown. The particle nature of 
the dark matter can be probed if it is detected either by direct detection process or by indirect 
detection. In the former the energy of recoil of a detector nucleus is to be measured if a dark 
matter indeed scatters off such a nucleus or nucleon. On the other hand the dark matter can be 
gravitationally trapped inside massive celestial objects such as sun, in regions of galactic centre 
etc. These trapped dark matter particles eventually annihilate to produce fermion-antifermion 
pairs or 7s. Study of dark matter through the indirect detection whereby such annihilation 
products are detected and analyzed, may reveal the nature of the dark matter. Recently much 
interest is generated by the evidence of monochromatic 130 GeV 7-rays from the direction of 
galactic centre obtained from the analysis [H [2] of observational data [3] from the Fermi-LAT 
satellite borne telescope. This phenomenon could not be explained due to other astrophysical 
events and the observance of 7-ray, is strongly suggestive of being produced by the annihilation 
of dark matter at the Galactic centre region. 

In the present work we propose a particle physics model for dark matter that can explain 
the 130 GeV 7-line observed by Fermi-LAT. Our model is in fact a two-component dark matter 
model in which a real scalar singlet and an inert doublet are added to the Standard Model (SM) 
of particle physics. There are previous works where either the real scalar singlet model or the 
inert doublet model has been discussed as one component dark matter model. But, as it will be 
revealed later that any one of these two models fail to explain individually (as one component 
dark matter) the observed 130 GeV 7-line from the galactic centre if it is produced due to dark 
matter annihilation. 

In earlier works such as Ref . 0j , the authors showed that the scalar dark matter can annihilate 
into 77 final state with the help of additional charged scalars in a model independent way. It 
would yield 130 GeV 7s with the required annihilation cross section ((<tv 77 ) ~ 10 _27 cm 3 /s given 
by the analysis pQ of the Fermi-LAT data mentioned above). Different other possibilities involving 
new particles originated from different extensions of SM have been investigated jS] E] to explain 
the 130 GeV 7-line through DM (single candidate) annihilation. However most of these endeav- 
ours are restricted in obtaining (crv 77 ) ~ 10~ 27 cm 3 /s without going into detailed discussions on 
relic density calculations, calculation of scattering cross sections relevant for direct detection and 
their comparisons with experimental (direct detection experiments) or observational (WMAP 
[7]) results. 

One of the simplest choice to accommodate a dark matter is to extend the SM with a gauge 
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singlet real scalar field S (we will call it real scalar singlet dark matter model, RSDM, from now 
on) [H 13 [10], which couples to SM Higgs (h). The use of a Z 2 symmetry ensures the stability 
of the dark matter candidate. In ref.[JJj] it is shown that this type of model fails to explain the 
130 GeV gamma-ray line from the Galactic centre. The reason of this failure is related with the 
relatively small annihilation cross section of two S fields into two 7s considering the fact that 
S should contribute to the correct amount of relic density as predicted by WMAP data [7]. In 
[9], it was shown that with SM Higgs much heavier than 125 GeV could in principle lead to 
(crv 77 ) ~ 10 _27 cm 3 /s with mg ~ 130 GeV, when the constraint on the S field to produce right 
amount of dark matter abundance is relaxed. However once this constraint is applied, the (crv 77 ) 
becomes few orders of magnitude less than required. In addition to these findings, if we employ 
the recent constraint on the Higgs mass from LHC experiment [11 J on the SM Higgs, then we 
infer that strictly within RSDM picture, we can not accommodate both the relic density as well 
as the 130 GeV 7-line from DM annihilation. 

Another well motivated dark matter model is the inert doublet model (IDM) [121 [TBI [TT] 
which requires an extension of the SM by a scalar Higgs (inert) doublet <3> having a Z 2 . In 
ref.[J_2], it was shown that there exists an allowed region (consistent with the WMAP results of 
relic density) in IDM for dark mattei§| 0° having mass in the range between 80 GeV and 160 
GeV provided the mass of the dark matter candidate (m^o) is less than the mass of SM Higgs 
(and top quark), m^o < mh,t- This condition was imposed to reduce the (av) total (to get rid 
of the contributions like 0°</>° — > hh (and/or tt)) since otherwise the relic density would be 
smalj§. Then they found that due to accidental cancellations of different Feynman diagrams for 
annihilation into gauge bosons, (crv) to tai can indeed be the right amount for a judicious choice of 
parameter space involved in the model. 

We are trying to find a resolution where a right amount of dark matter relic density could 
be obtained as well as an explanation of the 130 GeV gamma-ray line of Fermi-LAT can also be 
probed through DM annihilation. Now with the consideration that the SM Higgs boson has a 
mass rrih ~ 125 GeV [TTJ and mass of the DM candidate 130 GeV (in order to explain the 130 
GeV gamma-ray line from Fermi-LAT, the mass of the DM should be in this range), the above 
mentioned condition m^o < irti lt related to the IDM is evaded thereby the channel of annihilation, 
<p (j) — > hh, opens up. Hence (cv) total will increase and the final relic density in this sort of model 
can accommodate only (10-30)% of the observed DM relic density [13]. However this result has 
an interesting consequence. We can compensate this deficit of the DM relic density by another 
candidate of DM while 0° explains the 130 GeV gamma-line through its annihilation. 

Keeping in mind the above mentioned scenarios (particularly the RSDM and IDM), we pro- 

5 (f>° is the neutral component of the extra Higgs doublet. 
6 since relic density Slh 2 oc l/(crv). 
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pose that the dark matter can actually be composed of two fields, namely a scalar singlet (S) 
and an inert doublet ($). The additional feature of this model would be that these two compo- 
nents possess an interaction between them via a term like ($^$)5'S'. The 130 GeV 7-line will be 
produced by the annihilation of the component </>° while the role of other component S, besides 
contributing to the overall relic density is to increase 0° contribution to the combined relic den- 
sity through ($t$)S'S'. A model with a scalar singlet and a doublet was presented in [15] in the 
context of GUT models. A multi-component dark matter was proposed in [16] earlier where an 
additional fermion singlet and a scalar singlet were introduced, though not in the context of the 
possibility of having 130 GeV gamma-ray from DM annihilation. We have imposed a discrete 
Z2 x Z' 2 , under which S and $ transform non-trivially. 

The calculation of the flux for 130 GeV 7-ray from galactic centre region also requires the 
knowledge of dark matter density in the region of the galactic centre. In the absence of a unique 
density profile in literature we consider in the present work two dark matter density profiles 
namely NFW (Navarro-Frenk- White) profile [T7] and Einasto profile [18] and compute the flux 
using the cross section (crv^o ( ^o_ s>77 ) calculated in this work from our model. 

The paper is arranged as follows. In Section 2 we describe the structure of our two component 
dark matter model. Section 3 discusses the calculations of the relic densities of each dark matter 
components and hence the combined relic density by simultaneously solving the two Boltzmann's 
equations for the two components. In Section 4 we discuss how to constrain the parameter space 
of the model by comparing the relic density results with WMAP and the scattering cross section 
results with the direct detection experiments data. Section 5 gives the cross section calculations 
for 130 GeV 7-ray from annihilation and hence the calculations of 7-ray flux using different 
dark matter halo models. Finally in Section 6 we present some discussions and conclusions. 

2 The Two Component Dark Matter Model 

In the present work, the DM sector is composed of a real gauge singlet scalar field, S and an 
extra (in addition to the usual Higgs doublet, H) scalar doublet field, $ (doublet under SU(2) L ). 
An exact (unbroken) discrete symmetry Z 2 X Z' 2 is imposed under which all the SM particles 
are even, i.e. having Z2 x Z' 2 charge as (1,1) and for $, S the Z2 x Z' 2 charges are (1,-1) and 
(-1,1) respectively. Thereby both S and $ are fermiophobic and do not develop any vacuum 
expectation value (VEV) and hence inert. The construction therefore ensures the stability of 
both S and $ in this two component dark matter scenario. The scalar doublet $ can be written 
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as 

$ = ( /hao ) . (1) 

The Lagrangian of the construction can be read as, C = Am+Adm, where Am is the Standard 
Model Lagrangian and £dm stands for the DM sector consistent with all the symmetries. £dm 
is then given by, 

Adm = £rsdm + Adm + Ant, (2) 
where £rsdm and Adm re f er to the individual Lagrangian of real singlet scalar and inert doublet 
dark matter respectively and Ant is the additional part representing the interaction between 
the two components of DM. So the most general form of them, consistent with the SM gauge 
group as well as the discrete symmetry imposed, is as follows, 

£rsdm = ^^S-|U 2 -^S 4 -A 6 (HtH)SS, (3) 

Adm = (D„<f>)\D^) -P2(&®) 2 -Xi{tfm&®) - 

A 2 ( < & t H)(H t $) - A 3 [($ t H) 2 + h.c] , (4) 

Ant = -\ 5 (&$)SS. (5) 

The SM Higgs Lagrangian is included in the Am- In this model we have five new particles, 
two charged scalar (0 ± ) and three neutral scalar particles (0°, S, A ). Due to the stability and 
electrical charge neutrality we consider 0° and S as two viable components of dark matter in 
this model. 

After spontaneous breaking of the SM gauge symmetry, the masses of these new particles and 
Higgs are given by, 

m l± = i4 + ^ A ^ 2 ' ( 6 ) 

2 2,2 
m 4P = 1*2 + aV , ( 7 ) 

m 2 A o = ti 2 2 + (3v 2 , (8) 
m 2 s = k\ + \ 6 v 2 , (9) 
m 2 h = 2 Pl v 2 , (10) 

where v (= 246 GeV) is the Higgs VEV and p\ (= \j ^) is the coefficient of the quartic coupling 
of the SM Higgs (part of Am here) potential. The parameters a and (3 are defined in terms of 

A's as 

ol = ^(A! + A 2 + 2A 3 ), 

(3 = ^(A! + A 2 -2A 3 ). (11) 
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As is evident from the above set up, the model involves 10 parameters in total, specifically rrih, 
m^o, rriAo, m^+, ms, a, A5, A6, K2, Pi- Now fixing the Higgs mass mu at 125 GeV, we have 
altogether 9 parameters, which are further restricted from theoretical bounds as well as from 
experimental results as discussed below. 

• Vacuum Stability - The Lagrangian of this model (Eq. (j2J)) must be bounded from below. 
This condition will be satisfied if 

Pi, P2, k 2 > , (12) 
a, (3 > -VP1P2 , (13) 

Ai > -i^fpm, (14) 

A 6 > -2^/pJi^, (15) 
A 5 > 2 V /^. (16) 



• Zero VEV of $ and S - Ground state of the Lagrangian Eq. (T2J) must preserve Z 2 x Z' 2 
symmetry for stability of the dark matter candidates, this leads to the condition that the 
VEV of both $, S is zero. 

• Perturbativity- In order to be within the perturbative limit, the model parameters cannot 
be too large. This can be ensured provided 

|model parameters) < 4n. (17) 



• Neutral Scalar Mass - The LEP [19] measurement of the Z boson decay width leads to 
the condition 

m^o + m^o > m z . (18) 

• WMAP Limit - The Combined relic density of the dark matter components must satisfy 
the WMAP limit [7], 

0.1053 < tt DM h 2 < 0.1165 at 68% C.L. (19) 

for the dark matter in the Universe. This condition will further constrain the parameter 
space of this model discussed above. 

• Direct detection limits of dark matter - The results of the ongoing experiments for 
direct detection of dark matter also impose additional limits on the relevant parameters of 
the present two component dark matter model. Being inert, the elastic scattering between 



6 



N 



N 



Figure 1: Feynman diagram for the elastic scattering between dark matter particle £ (i,e. S, 4>°) 
and nucleon N of the detector material via Higgs exchange. 

the dark matter candidate £ (here 0° and/or S) and nucleons (N) can take place only with 
the SM Higgs exchange as in Fig [TJ The relevant term in the Lagrangian which describes 
the interaction between £ and Higgs is given by 

£ = -k£ 2 h, (20) 

where the coupling k = av (A 6 t>) for £ = 0° (S). The spin independent scattering cross 
section for this process, £iV — > £iV, is given by[20j, 

- SI = ^f-V^. (2D 



7r \v J m| m h 



where and are the masses of the DM candidate and nucleons respectively, /i 



is the reduced mass and f represents the strength of the effective interaction which 
depends upon the number of heavy quarks involved [21J. In this work, we have considered 



the value of / = 0.3 

An upper limit on <r SI for a particular mass of dark matter particle would automatically 
sets an upper bound on the absolute value of the couplings a and A 6 . Since our goal is 
to explain the 130 GeV 7-ray line from annihilation of one of the DM candidates in the 
present two component dark matter model, that particular component should have a mass 
near 130 GeV. The upper limit on cr SI for a 130 GeV DM candidate as obtained from 
XENON 100 (2012) [23] data, is ~ 3 x 10" 45 cm 2 at 90% C.L., which in turn specifies the 
range of a and A 6 as 

H , |A 6 | < 0.038, (22) 
The similar ranges for a, \$ from XENON 100 (2011) [21] data are obtained as 

\a\ , |A 6 | < 0.076. (23) 
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The result obtained from XENON 100 (2011) data was less constrained than 2012 data. 
This will be discussed in Section 4 while constraining the parameter space. Upper limits 
set by XENON 100 (2012) on <r SI with respect to dark matter mass lie below the similar 
limits obtained from other recent ongoing dark matter direct detection experiments such 
as CDMS-II [25J, EDELWEISS-II [26] etc. Therefore the limit given in Eq. ([22]) for a 130 
GeV dark matter particle must also satisfy the limits (from (cr SI — m x ) exclusion plots) 
obtained from other direct detection experiments . 



3 Combined Relic Density Calculation for the two dark 
matter Candidates 0° and S 

In the present two component dark matter model, the total relic density of the dark matter in the 
universe will have contributions from both the components S and 0°. While both are annihilating 
into the SM particles, the heavier component, S can annihilate into the lighter component, 0° too. 
In order to obtain the correct combined relic density we have to solve Boltzmann's equation for 
each components simultaneously. The coupled Boltzmann's equations [27] to study the evolution 
of the number densities of the two dark matter candidates (ns and n^o) are given by, 

^ + 3n s H = -(ayss^xx) ( n s ~ ( n s) 2 ) ~ (<™ss^4P4fl) (ri 2 s - j > ( 24 ) 



dn^o 

— + 3n,oH 



( \ / (n eq ) 2 \ 

-(ctv^o^) (nj - (n e fi 2 ) + (o-vss-^o) I n% - 7^fp n J° 1 • ( 25 ) 



Here n e ^ and n e g are the equilibrium values of n^o and ns respectively, H is the Hubble's con- 
stant, x represents any SM particle such as leptons, quarks, gauge bosons, Higgs boson. The 
annihilation of SS into is included through the annihilation cross section ((rvgg^o^o), the 
expression of which in our scenario is discussed later in this section explicitly along with the total 
annihilation cross section. 

Introducing two dimensionless variables Yi — — and Xi = ^ with i = S,<f)°, where s and 
T are the entropy density and temperature of the universe respectively, Eqs. 02U can be 
written as 

dY s f45G\-tm s ^f, , fv3 ™, 2 , , , . {Yl q f 



dxf 



(26) 
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dY^o 



45G 



7T 



9* ( (<™<t>°4> Q ->xx) ~ ( Y £o) 2 ) ~ (vvss^<p°) 



(Ys Q ) 2 



(27) 



Here G is the Gravitation constant, and g* is denned as, 

Mr) / i din(Mr)) \ 

^ = 7swl + 3^mttJ • (28) 

with g e s{T) and h e s(T) are the effective degrees of freedom related to the energy and entropy 
densities through p = # eff (T)|^T 4 , s = /i c ff(T)^T 3 . 

Once we get the values of Y^o and Ys at the present temperature T after solving the coupled 
Eqs. (f26 | |27|) . we will be able to calculate the individual contributions and Qs from [28] . 



a/i 2 = 2.755 x 10 s 



GeVJ 



) Yi(T 



(29) 



(using the present values of s and h). In this work we have solved the coupled Boltzmann 



equations numerically to get the values of Yi(T ) 
density of the universe can be obtained through 

Q r h 2 = VLjjah?' - 



After having these estimates, the total relic 



(30) 



As we mentioned before, (avss-+ X x) ^ n Eq. (l2"4"|) represents the total annihilation cross section 
of two S particles into SM particles such as leptons and quarks (//), gauge bosons (H^ + H^~, ZZ) 
and Higgs boson (h). The Feynman diagrams for all the processes are shown in Fig. [2] The 
expressions of annihilation cross sections of two S particles for these final states are given below 
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Figure 2: Lowest order Feynman diagrams of two S annihilate into a pair of fermion and anti- 
fermion, W + W~ , ZZ and Higgs. 



In writing the expression for the cross sections, we have introduced the notation g xyz (or g xyzp ) 
which are related to the coupling constants of the corresponding interaction Lagrangian through 
£ — 9xyzXY Z (or C = g xyzp XYZP). Here X,Y,Z,P are the fields involved in a particular 
process. All the gxyz and gxyzp are listed in Table [U The masses for fermion /, W boson, Z 
boson and Higgs boson are denoted by m/, mw, mz, rrih respectively and Th represents the Higgs 
decay width. 

Similarly in Eq. (l25p. (crv^o^^) represents the total annihilation cross section of two <p° 
particles into SM particles. The Feynman diagrams of individual processes are shown in Fig. |3] 
and the expressions of annihilation cross sections for these channels are given below. 
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interactions involving XY Z{F) 
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Table 1: Interactions and the corresponding couplings. 
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Figure 3: Lowest order Feynman diagrams of two 0° annihilate into a pair of fermion and anti- 
fermion, W + W~ , ZZ and Higgs. 
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The other cross section involved in the coupled Boltzmann equations is (avss->< 
to the presence of the interaction between the two dark matter components S and 
expression for this cross section is given by 
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The heavier component of DM is chosen to be S in our casqj. Since only will contribute to 
the 130 GeV gamma ray production (the annihilation cross section ((rvss-+ry) for the singlet S 
is not enough to produce Fermi-LAT observed 130 GeV 7— ray flux as we argued before), it's 
mass should be m^o ~ 130 GeV. We will show later that <p° with its mass ~ 130 GeV can not 
accommodate the total relic density consistent with the WMAP data. So the rest of the relic 
density should be provided by the other component S. And we find that it would be a good choice 
to make mg > m^o. The involvement of this interaction between two dark matter components 
is a salient feature of our analysis. It would enhance the number density of 0° through the 
annihilation of SS during evolution. The Feynman diagrams for this process (SS — > 4>°<f) ) are 
given in Fig. HJ 



5- _ -c 



^ If 

s. ✓ 

X 



Figure 4: Feynman diagrams for the annihilation channel SS — > 



In the left panel of Fig. [51 we have shown the variation of (av^o^o^^) with parameter a (Eq. 
(PH) f o r different values of m A o, namely m A o = 600 GeV, 500 GeV, 400 GeV, 350 GeV, 300 
GeV. In the right panel, the variation of (crvss^xx) "with parameter A 6 (Eq. Q) is shown. We 
have chosen m^o = = 130 GeV, a = -0.037 and mj = 130.5 GeV for drawing these plots. 
The reason behind this choice of parameters should be cleared in the next section. 



4 Calculation of Model Parameters 

In this section we describe the procedure adopted in this work in order to estimate the values 
of the parameters involved in our proposed model. At the very outset, we summarise few basic 
requirements. In the present two component dark matter model, only the component <j)° can 
account for the observed 130 GeV 7-line through the annihilation channel — y 77 details of 
which is discussed in the next section. This is due to the fact that the other component S being 
a scalar singlet cannot produce sufficient annihilation cross section for the channel SS —y 77. 

If ms < rn^o, (jp component would be annihilated into S (via the process 4>°(j> — > SS). This would reduce 
the contribution of <fi° to the combined relic density and hence the 7— ray flux originated from </>°^>° annihilation 
would be suppressed. 
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Figure 5: Left panel - Variation of (ctv^o^^) with a for different values of m^o, Right panel 
Variation of (crvss-^xx) with A 6 for m s = 130.5 GeV. 



So we conclude that the mass of the 0° component needs to be ~ 130 GeV. We also find that a 
choice for the mass of the charged scalars, <f> (involved in the loop of the 0°0° — > 77 process) 
close to enhances the annihilation cross section (crv^o^^) (see Eq. (f4"oT) ) . Therefore . in order 
to maximise this contribution (required to achieve the cross section in the right ball park without 
taking a very high value for the relevant parameter in the model), we consider m^o ~ m^,± = 
130 GeV throughout the present discussion. Note that such a choice is consistent with the LEP 
bound [301. With this consideration, the parameter a is reduced to Ax/2. 

Of course 0° with m^o = 130 GeV cannot individually account for the WMAP results on 
relic density. In our scenario, this deficit will be compensated by the contribution of the other 
component S. However we need to maximise the contribution of 0° to the combined relic density 
so as to keep the flux of the observed gamma ray (130 GeV) from the Galactic centre at an 
adequate level (see Eqs. (ll2l[B^]) ). For example if 0° contributes to 60% compared to a case where 
it contributes only 30% to the total relic density then the flux for the gamma ray originated from 
DM annihilation will also be proportionately higher compared to the latter case. Apart from the 
WMAP data, we also use the limits obtained from the dark matter direct detection experiments. 
Needless to mention that in doing so the conditions obtained in Eqs. (U2H16t W7\ [TH]) are always 
satisfied. 

In order to find a choice of parameter space for which the contribution from 0° toward relic 
density can be maximised, we calculate the ratio ( Q° h 2 ) ■ The relic densities Vt^oh 2 and Vt c h 2 
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we mean that 



are computed using Eqs. fl26| 1271 129|) . In the left panel of Fig. [6] we plot ^7^5- 
function of the parameter a for different values of m^o. Here by ( q°^2 ) j 

V c / Maximum 

we have chosen only the largest possible value of the ratio corresponding to a particular a while 
other parameters are scanned over their entire range (—0.1 < A6 < 0.1, 0.1 < A5 < 1.0). The 
limits of a (Eqs. (|22l[23])) obtained from XENON 100 (2012) data and XENON 100 (2011) data 
are superimposed on this figure (left panel of Fig. [6])) to find out the allowed region for the ratio, 
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Figure 6: Left panel - Variations of 
Right panel - Variations of 



n c h 2 



n c h 2 



Maximum 



with parameter a for different values of m^o, 

Maximum 

with mjsfi for different values of parameter a. 



In the left panel of Fig. Othe slategray colour band represent the limits on a (\a\ < 0.076) 
obtained from XENON 100 (2011) data while the limits on a (Eq. [22D from XENON 100 (2012) 
data are depicted in the figure with the turquoise colour band. We want to find a suitable value 
of a (arguably from the region allowed by XENON 100 (2012), as then the other bounds will also 
be satisfied automatically) for which we get a maximum contribution to the ratio ' 



Q c h 2 



Note 



we show the variations of 



that this ratio depends on the choice of m^o. In order to have more insight on this dependence, 

^° h 2 ) with mfio for different values of a (which are allowed 

c / Maximum 

by XENON 100 (2012) data) in the right panel of Fig. [6j It is evident the contribution of (f>° 
increases as m^o increases. This can also be noted from the same figure that for a fixed value 
of a, the increase of Q^oh 2 with m^o is stripper for lower values of m^o. However we cannot 
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choose a value of m^o which is arbitrarily large since our choice m^o = m^± = 130 GeV imposes 
a relation between parameters A2 and A3. For example, to achieve m^o = 500 GeV, we find 
A 2 = 3.852, A 3 = —1.926 when a = -0.037 is considered (fixing A 2 and A 3 will be discussed 
later in the context of Fig. [7]). Taking into consideration all the above discussions we can have 
the contribution of 0° to the total relic density to be at most ^ 62%. A choice of m^o = 900 
GeV would require a value of A 2 ~ 47r, which poses a threat to the perturbativity. Therefore 
in the present work we consider m^o < 500 GeV and the corresponding parameters as obtained 
following the discussions above. 




4 5 6 7 
Am (m s -m ») in GeV 



0.6 



0.5 







a 0.4 



0.3 





























































































































































































































m s = 130.5 GeV 

m s = 131.0 GeV 

m s = 132.0 GeV 






















m s = 133.0C 
m s = 135.0C 


eV 

eV 

eV 












m, = 140.0 C 



150 200 



250 



300 350 



400 



450 500 



Figure 7: Left panel - Variations of 
m^o, Right panel - Variations of 



n c h 2 
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with Am(ms — m^o) for different values of 

Maximum 

with rriAo for different values of Am(ms — m^o). 



As it is now evident that the inert doublet component alone with m^o = 130 GeV cannot 
account for the total dark matter content of the universe and one needs to add the relic density 
of the singlet component 5* for producing the WMAP satisfied total relic density. In fact this 
is one of the motivations for choosing this two component (inert doublet + scalar singlet dark 
matter model). In order to choose a suitable mass for S in the present two component dark 
matter model, we define a quantity Am, (= mj — m^o) remembering that ms should be heavier 



than m^o as discussed before. We then study the variations of 



n c h 2 



with Am for 



Maximum 

different values of m^o which is shown in the left panel of Fig. [7J In the right panel of Fig. [7] 
we show the variations of the fraction ( -^5 ) with m° A (< 500 GeV) for different values 

V c / Maximum 

of ms- The plots in both the panels are obtained with a = —0.037 (which satisfy both XENON 
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data). Comparing both the panels of Fig. [7] one concludes that the contribution of </>° to the 
combined relic density increases as the mass splitting between the two components of the dark 
matter decreases and the ratio ( ^° fe2 ) is ~ 62% for Am = 0.5 GeV and m A o = 500 

V c / Maximum 

GeV. 




2e-26 4e-26 6e-26 8e-26 1e-25 1.26-25 1.4e-25 1.6e-25 1.8e-25 1e-30 1e-29 le-28 1e-27 1e-26 1e-25 

^smV* (cm3/s) <0V sa + xP (cm3/s) 

Figure 8: Left panel - Variations of the contributions of S (Qsh 2 ) and 0° (Q^oh 2 ) to the combined 
relic density (Q c h 2 ) with the annihilation cross section (<JVss^4> Q <f>°), Right panel - Variations of 
the contributions of <p° (Q^oh 2 ) to the combined relic density (Q c h 2 ) with the annihilation cross 
section (av ss ^ xx ). 

We also calculate the variations of (^fj^J and (^jf^j with the cross section (crygs^o^o) 
and the results are plotted in the left panel of Fig. [HJ In all these of calculations, we have 
used our standard set of parameters, m^o = = 130 GeV, mj = 130.5 GeV, a = -0.037, 
m^o = 500 GeV (as mentioned earlier). From this plot, it is clear that initially when ((tvss^</> ^°) 
is very small (or nearly zero) the contribution of 4>° is only ~ 25% of the combined relic density. 
This is because for small values of (crv.55_j.0O0o), Eqs. fl2"ET |2T|) effectively become two decoupled 
equations that represent the Boltzmaan's equations for RSDM and IDM respectively. Under 
such circumstance, the calculations of the individual contributions for S and 0° are pursued . 

Note that in this decoupled scenario, individual relic density contribution would be inversely 
proportional to the corresponding annihilation cross section. Therefore we need to have an 
estimate for (av^o^o^^) and (<jv S s^ X x)- ^ n or der to understand this in more detail, we refer to 
Fig. [5] From the left panel of Fig. [51 we see that the value of (av^o^o^^) is ~ 8.5 x 10 -26 cm 3 /s 
for rriAo = 500 GeV, a = —0.037, m^o = = 130 GeV. This value of (crv^o^^) is nearly 
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4 times larger than what is required to get WMAP satisfied relic density for dark matter mass 
of 130 GeV, thereby contributing only ~ 25% to the total relic density (Eq. (j!9p ) as seen from 
left panel of Fig. [SJ So we would expect that the rest ~ 75% contribution should come from the 
singlet scalar component S in the present model. This is indeed possible if we consider the right 
panel of Fig. |5j This corresponds to a case where practically there is no interactions between S 
and 0°, i.e. A5 ~ 0. As the interaction between S and 0° becomes increasingly stronger (i.e. |As| 
starts to have nonzero value), more and more S particles annihilate to produce 0° particles and 
contribution of 0° to the combined relic density will be boosted. It reaches a maximum which 
comes out to be 62% in our present analysis for (crvss->4>°(f> ) ~ 1-8 x lCT 25 cm 3 /s as seen from 
the left panel of Fig. |HJ 

The right panel of Fig. shows the variation of \^=fjjr\ with (crvss-+xx}- ^ indicates that 
contribution of 0° decreases as the value of (crvss^xx) increases. This is because of the increment 
of (cvss^xx) signifies large number S particles annihilate into SM particles and consequently 
less number of particles are available for annihilation of S to produce 0° in the final state. As 
a result the relic density of 0° decreases. From this plot, we conclude that in order to obtain 
the relic density contribution of 0° close to 62% or above the value of (<rv ss ^ xx ) should be 
^ 1.74 x 10~ 28 cm 3 /s. Note that (avss-*xx) involves the parameter Xq and A5, Xq both of which 
are involved in ((rvss-^^o). Therefore this limit (on (avss->xx)) a l° n g with the other one we 
just discussed above, (avss-+<j>°<]> ) ~ 1-8 x lCr 25 cm 3 /s, set up a range of allowed region of A5 and 
A6 if we restrict ourselves with (^fljr) ~ 62%. The allowed range of values of parameters A5, A6 
for the present case (m^o = 500 GeV) is shown in Fig. [9j In Table [2] we furnish the values of 
the parameters A5 and Xq along with other model parameters. 

0.3873 1 , , , , , , , 1 




Figure 9: Allowed range of the parameters As,A6 in order to obtain (o~vss^ X x) — x 
10 _28 cm 3 /s and (avss-*j>°4> ) ~ 1-8 x 10 _25 cm 3 /s such that the contribution of 0° component is 
maximum to the combined relic density for m^o = = 130 GeV, Am = 0.5 GeV, a = —0.037, 
m A o = 500 GeV. 
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Mass of 
A (m° A ) 
(GeV) 


Contribution of 
(f)° in the combined 
Relic Density 


^2 

(GeV) 


M 


^5 


I \ I 

l A 6| 


^2 


^3 


300 


51.7% 


138.344 


-0.074 


0.2698-0.2703 


< 2.6 x 10~ 3 


1.208 


-0.604 


350 


55.3% 


138.344 


-0.074 


0.3014-0.3020 


< 2.8 x 10~ 3 


1.745 


-0.872 


400 


58.1% 


138.344 


-0.074 


0.3322-0.3327 


< 2.8 x 10~ 3 


2.365 


-1.182 


500 


62.1% 


138.344 


-0.074 


0.3867-0.3873 


< 3.2 x 10~ 3 


3.852 


-1.926 



Table 2: Relic density contribution of <f)° for different values of m° A and the corresponding values 
of model parameters for a = —0.037, Am = 0.5 GeV (ms = 130.5 GeV, m^o = 130.0 GeV). 

5 130 GeV 7-ray line from Dark Matter annihilation 

In this section our endeavour will be to explain the recently observed 130 GeV 7-line from 
the Galactic centre originated from dark matter annihilation in the framework of the present 
two component dark matter model. In order to produce such a 130 GeV 7-line, the required 
DM annihilation cross section into two photons should be ~ 10~ 27 cm 3 /s as predicted from the 
analysis [H [2] of Fermi-LAT data [3] . In the context of the present two component dark matter 
model, only the <f)° component having mass 130 GeV can contribute to the production of 130 
GeV 7— ray line. The annihilation of </>°0° into 77 can take place only via charged scalar ^ and 
W ± loops. It can indeed produce the required cross section ~ 10 27 cm 3 /s. The cross sections 
for other annihilation channels that can produce 77 (e.g. via Higgs) for both the components 0° 
and S are orders of magnitude less than this value [TU] . 

We calculate the 7-ray flux due to annihilation of 0° in the "central region" of our Milky way 
galaxy. The lowest order Feynman diagrams for the process — > 77 via ± and W ± loops 
are shown in Fig. [T0J The expression of differential 7— ray flux due to dark matter annihilation 
in galactic halo is given by 



(i$ 7 1 (av^o^o^yy) dNy 2 



TqPqJ , (40) 



where r = 8.5 kpc is the distance of the sun from the Galactic centre and p Q = 0.4 GeV/cm 3 
is the the local dark matter density at the solar neighbourhood. The quantity J in the above is 
given by 

J = I dl I db cos 6 J{l,b) , (41) 
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Figure 10: Lowest order Feynman diagrams for the process 



_> 77 



with 



and 



J(l,b) 



ds ( p(r) 



i.o.s r e V Pe 



(42) 



AQ = A dl db cos b 



(43) 

Here / and b in Eqs. flJTJ S3]) are galactic longitude and latitude respectively. In the above, 
the integration is performed around a radius of 3° (galactic central region) around a centre with 
coordinates (I, b) = (—1° , —0.7°) [2J. In Eq. (j42j) r and s are related by, 

r = (s 2 + Tq — 2sr cos I cos 6) ^ 2 . 



The expression of the energy spectrum of 7, denoted by ^p- is given by, 



dN 1 
dK 



25{E - E, t 



(44) 



(45) 



We have performed /, b integration (in Eqs. (141 [ 143]) ) over the "central region" of our galaxy and 
the s integration (in Eq. |4"2"]) is along the line of sight (l.o.s). 

We explicitly calculate all the Feynman diagrams given in Fig. [10] and obtain the expression 

for (av^o^o^yy) as 



2 2 
at m^o 

32tt 3 



p°(f> <j)+<j)- r <j) 



Fj,+ 



i,°s w+w- r w 



m 



4>+ 



m 



(46) 
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Where 



with t 



4 m 7 



V = r[l-r/(r)] 
F w = 2 + 3r + 2r(2-r)/(r). 

!) ± ,iy ± , and 
sin - 

1 (V^7^)I < for r > 1 



(47) 



sin 

1 
4 



1 + V T^7 . 

in , — i7T 



(4£ 



for r < 1 . 



1 - 

Since the couplings gy^o ^+^,- = — p 2 and gy>0ow+w- = ^f 1 (Tabled]), Eq. fT46|) can be written as 



((TV $0 ^o^. 7J ) 



a 2 m| 
32tt 3 



P2F0+ F w 
2 + „,2 



0+ 



(49) 



In the present work, 7-ray flux is calculated for two different dark matter halo profiles namely 
the NFW profile [17] and the Einasto profile [18]. These halo profiles give the functional depen- 
dence of dark matter density p(r) with r. The expression of p{r) for the NFW profile is given 
by, 



PnfwIH 
and for the Einasto profile 

PEinasto(^) = P S exp 



Ps 



(r/r s ) (1 + r/r s y 



(50) 



r 2 




1 


1 7 







(51) 



where r s in Eqs. ( l5"Ull5"l"j) is taken to be 20 kpc and 7 = 0.17 in Eq. (|5l|) . In the above the 
value of the normalisation constant p s is determined by demanding that at r = r , the density 
p(r) = pq. We have seen earlier that in the present model of two component dark matter, the 



n c h 2 



0.62 ) of the total dark matter relic 



inert doublet component 0° contributes to ~ 62% 

density. Therefore in calculating the 7— ray flux from the process cfPcfP — > 77 we compute p s by 
taking 

O.n/l 2 

(52) 



Pq X 



tt r h 2 



and demanding that for the dark matter component <ft°, p(r) = p' Q at r = r with p & = 0.4 
GeV/cm 3 [32]. 
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Mass of 

A 
(Lrev ) 




P'e 

(GeV/cm 3 ) 


Ps 

(GeV/cm 3 ) 


(<7V«^ -+7<y) 

(cm 3 /s) 


P2 


500 


0.621 


0.248 


0.214 


i.i83i8:^ x io- 27 


^- '^-0.64 


400 


0.581 


0.232 


0.201 


1.3411°;^ x 10- 27 


c oo +0 - 67 

d - uu -0.68 


350 


0.553 


0.221 


0.191 


1.489+°;^ x 10- 27 


r 90+O.7O 
d - zo -0.74 


300 


0.517 


0.207 


0.178 


1.710j£f$ x 10- 27 


e cc:+0.75 
J.Jd_ 76 



Table 3: Results for the NFW Profile. 

From the left panel of Fig. 3 of Ref. [2] the best fit value of the 7— ray flux (in terms of 
E 2 <& (GeV cm -2 s" 1 sr" 1 )) observed by Fermi-LAT from the central signal region of the Galaxy 
can be read as E 2 <& = 5.6 x 10~ 5 GeV cm -2 s _1 sr _1 with 95% C.L. error band that lies in the 
range 3.97 x 10 -5 < E 2 § < 7.51 x 10~ 5 (GeV cm" 2 s _1 sr _1 ). We use these values of the flux in 
Eq. PHI) and compute (av^o^o^^) for the best fit value of the flux as also the two extremities 
of its error band at 130 GeV. Note that the parameter P2 is yet to be determined. This can be 
estimated by calculating the cross section (av^o^o^^) (given by Eq. (j4"9l ) and hence the flux 
E 2 § and then comparing this flux with that given by the Fermi-LAT data. 

The results are furnished in Table [3] and Table H] for the NFW profile and the Einasto profile 
respectively. They are given for the chosen mass m^o of 500 GeV as also for three other values 
of rriAo namely 400, 350, 300 GeV for the purpose of demonstration. In both the Tables El H] the 
values of the cross sections obtained for the best fit value of the flux are given. The computed 
cross sections for the two extremities of the error band of the flux are shown by the subscripts 
and superscripts of the central values. The corresponding values of p 2 that are calculated using 
Eq. f l4"9]) are also shown in similar fashion. It is seen from both the Tables that although the 
calculated values for p 2 depend on the dark matter density profile that one chooses, they are 
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Mass of 

A 
(Lrev ) 




P'e 

(GeV/cm 3 ) 


Ps 

(GeV/cm 3 ) 


(cm 3 /s) 


P2 


500 


0.621 


0.248 


0.051 


0.608±g;?^ x IO" 27 


o ^7+0.47 
°- ' -0.50 


400 


0.581 


0.232 


0.047 


0.712l°JS x 10- 27 


q OO+0.51 
°- oz -0.52 


350 


0.553 


0.221 


0.045 


0.7801S >< !0" 27 


q 07+O.53 
°- JI -0.54 


300 


0.517 


0.207 


0.042 


0.896^ x 10- 27 


A O1+0.56 
^• zi -0.58 



Table 4: Results for the Einasto Profile. 

within the perturbative limit and the corresponding cross sections (o - v^o^o_ >77 ) are also within 
the desired limits of ~ 10 -27 cm 3 /s. 

Here we like to mention that we have checked the possibility that the continuum gamma- 
rays may overshoot the monochromatic gamma-ray line. In Ref. [6], the authors have shown 
that in order to distinguish the monochromatic gamma-line (from DM DM — > 77 channel) from 
the continuum gamma-ray spectrum (produced by the secondary photons originating from the 
annihilation products of dark matter e.g. gauge bosons, qq, //), the branching ratio for the 
channel DM DM — > 77 must be greater than 1% of total annihilation cross section (sum of 
annihilation cross sections for all possible channels). For the dark matter component 4>° in our 
model we find that this ratio is nearly 1/70 (i.e. > 10~ 2 ) as seen from the left panel of Fig. |5] 
and Table |3j 

6 Discussions and Conclusions 

In the present work we propose a dark matter model which contains two dark matter candidates. 
Such a two component dark matter model can be obtained by adding a scalar singlet S (singlet 



23 



under SM gauge group) and a doublet $ (doublet under SM gauge group) to the scalar sector 
of SM. We have introduced discrete symmetry Z2 x Z' 2 under which only S and $ transform 
non-trivially. Both the scalar singlet S and doublet $ do not produce any VEV. Consequently 
Z 2 x Z' 2 symmetry remains unbroken which ensure the stability of both the components (S, 4>°) 
of the dark matter in the present model. While the component 0° (neutral part of the doublet 
$) can produce the annihilation cross section required to obtain 130 GeV 7-line, the value of 
the corresponding cross section for the scalar singlet component S falls deficit by few orders of 
magnitude. However the component 4>° above, having a mass of 130 GeV cannot solely account 
for the relic density predicted by WMAP. This deficit in relic density is compensated by the 
scalar singlet component S such that the combined relic density (Q c h 2 ) for this two component 
dark matter model always lies within the range given by WMAP. Combined relic density is the 
sum of individual relic densities of both the components S and </>° which are obtained by solving 
the coupled Boltzmann's equations numerically. We have found that the contribution of the 
component 0° to Q c h 2 will be maximum (~ 62%) when we consider a = —0.037, | Xq\ < 3.2 x 10~ 3 
(both satisfy XENON 100 (2012) limit as well as limits from other dark matter direct detection 
experiments namely CDMS-II, EDELWEISS-II etc.), Am = 0.5 GeV, A 5 ~ 0.387 and m A o = 500 
GeV. Finally in the last section we have calculated the annihilation cross section (av^o^o^.^) for 
the channel cjPcfP — > 77 with the mass of <f)° ~ 130 GeV. Using the expression of this annihilation 
cross section ((crv^o^^)) we have computed the 7— ray flux of energy 130 GeV for two different 
dark matter halo profiles namely the NFW profile and the Einasto profile. The exact dark matter 
density at the galactic centre is unknown (e.g. Ref. [21] and references therein). This may 
produce an additional uncertainty in the flux calculation. Depending on the value of the dark 
matter density at the galactic centre, 130 GeV gamma-line may also be produced for a value of 
annihilation cross section lower than the specified value of ~ 10~ 27 cm 3 /s. Fermi Collaboration 
placed an upper limit [22] crv 77 < 1.4 x 10~ 27 cm 3 /s for 130 GeV dark matter with an NFW 
profile and <rv 77 < 1.0 x 10 -27 cm 3 /s with an Einasto profile. In the present work we indeed 
obtain av 77 in the range of this upper limit (~ 10~ 27 cm 3 /s). 
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